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Generation of few-cycle optical vortex pulses is challenging due to the large spectral bandwidths, as most 
vortex generation techniques are designed for monochromatic light. In this work, we use a spiral phase plate 
to generate few-cycle optical vortices from an ultrafast titaniurmsapphire oscillator, and characterize them in 
the spatiotemporal domain using a recently introduced technique based on spatially resolved Fourier transform 
spectrometry. The performance of this simple approach to the generation of optical vortices is analyzed 
from a wavelength dependent perspective, as well as in the spatiotemporal domain, allowing us to completely 
characterize ultrashort vortex pulses in space, frequency, and time. 
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1. Introduction 

An optical vortex is a field with a phase singularity, such 
that the phase is undefined in the center of the vortex 
and thus the intensity must be zero. In addition, the 
phase variation along a path enclosing the singularity is 
non-zero: 

i = ^/ v ^( r ) ' dr ' 0 ) 

where £, which must be an integer different from 0, is 
called the topological charge. In the simplest case, the 
phase front of the light twists helically like a corkscrew 
around the singularity, according to 0) = £0 , where 
p and 0 are the radial and angular components in cylin¬ 
drical coordinates. 

Phase singularities in propagating wave fields were for 
the first time observed by Nye and Berry in ultrasound 
fields in the 1970s [1]. Vortices in laser beams were later 
described [2] and it was recognized by Allen and co¬ 
workers in 1992 that helically phased beams carry orbital 
angular momentum (OAM) [3, 4]. The particular prop¬ 
erties of optical vortex beams have led to many applica¬ 
tions [5], such as imaging of faint objects around a bright 
one [6], overcoming the diffraction limit in stimulated 
emission depletion microscopy [7], and optical tweezers 
[ 8 ], 

Several methods exist today to convert a Gaussian 
laser beam into a vortex beam. The most commonly 


* Corresponding author: miguel.miranda@fysik.lth.se 


used are based on phase masks which are simply inserted 
into the beam in transmission [9], diffraction holograms 
[10], and refractivediffractive elements [11, 12]. Recent 
advances on spatial-light modulators enable the produc¬ 
tion of optical vortices with added flexibility, such as eas¬ 
ily changing target wavelength and topological charge. 

Recently, optical vortices have attracted interest in 
the ultrafast community [13-17]. Combining vortex 
pulses with high intensity ultrashort sources opens the 
way for strong-field physics with singular light beams 
[18, 19]. Generating ultrashort vortices is challenging 
since most techniques used to create monochromatic vor¬ 
tices have limitations when applied to broadband laser 
sources [20]. Even though polychromatic methods have 
been suggested and applied to ultrashort vortex gener¬ 
ation [13, 14, 21], these are usually more complex to 
implement, making monochromatic techniques still at¬ 
tractive. 

An easy way to create ultrashort vortex fields is to use 
a spiral phase plate, e.g. made of glass. Transmission 
through a material with thickness L leads to an optical 
phase of p(X) = (2ttc/ A)An(A)L, where A is the wave¬ 
length and c the speed of light in vacuum, and An the 
difference between the refractive index in the material 
and in vacuum. This phase varies with wavelength even 
in the absence of material dispersion (normal material 
dispersion makes the problem even worse). Therefore, 
the thickness jump of the spiral phase plate corresponds 
exactly to a phase equal to a multiple of 2i r only at the 
design wavelength of the phase plate. More generally 
the phase introduced by the glass plate varies according 
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Fig. 1. Simplified experimental setup. The beam from an ultrafast oscillator is compressed with a chirped mirror and fused 
silica wedge compressor. The tunable dispersion of the compressor is also used for the temporal characterization of the reference 
pulse with the d-scan technique. A dispersion-balanced interferometer, based on a MachZehnder design, is used to interfere 
a spatially-filtered reference wave with the generated vortex pulse on the chip of a CCD camera. The delay between the two 
fields is scanned and an interferogram is recorded for each pixel. 


to: 

v(p,e,\)=i(\)o. ( 2 ) 

It is therefore not clear what impact the chromaticity of 
the phase plate has on the quality of the (broadband) 
ultrashort vortices. 

Characterizing ultrashort vortex pulses is also chal¬ 
lenging: a few-cycle optical vortex is an inherently three- 
dimensional object in space and time, so its complete 
characterization requires spatiotemporal characteriza¬ 
tion techniques. Knowing the spatiotemporal evolution 
of the electric field of optical vortices is important for 
simulating nonlinear effects driven by high-intensity vor¬ 
tex pulses, such as HHG [22], as the effect is directly 
dependent on the electric field. 

In this work, we generate few-cycle optical vortex 
beams from an ultrafast titaniurmsapphire (Ti:Sapph) 
laser oscillator using a spiral phase plate and char¬ 
acterize them in the spatiotemporal domain, using a 
technique based on spatially-resolved Fourier transform 
spectrometry [23, 24]. This technique is well suited 
for this particular case where high spatial resolution is 
needed, since scanning is performed only in one dimen¬ 
sion (time). Our method allows for sub-cycle precision 
on the determination of relative phases across large ar¬ 
eas. 

2. Experimental Setup 

The experimental setup is schematically shown in Fig. 1. 
Three main parts can be distinguished: generation of 
the vortex pulses, spatiotemporal measurement, where 


the vortex pulses are characterized by interfering them 
with a reference, and characterization of the reference 
pulse. The ultrafast oscillator (Venteon Pulse-one) used 
as laser source for the experiments has a repetition rate 
of 80 MHz and around 3 nJ of energy per pulse with 
a spectrum supporting 5.0 fs pulses. Chirped mirrors 
(Venteon DCM7) provide negative dispersion, and fused 
silica wedges are used for dispersion fine-tuning. A 1 mm 
thick spiral phase plate (RPC Photonics) with topologi¬ 
cal charge of one at 800 nm, lithographically etched in a 
fused silica plate, is used in transmission to convert the 
laser output to a vortex beam. The dispersion of the 
plate can easily be compensated for by removing glass 
insertion from the wedges. Unlike approaches based on 
polymeric phase plates or spatial light modulators, the 
etched phase plate can stand high intensity and high 
average power and is well suited for applications in non¬ 
linear optics and strong-field physics. 

The spatiotemporal characterization is performed us¬ 
ing a technique based on spatially-resolved Fourier trans¬ 
form spectrometry [23]. Part of the original beam is 
split before the spiral phase plate, sent through a delay 
stage, and used to generate a homogeneous spherical ref¬ 
erence wave by focusing with an off-axis parabolic mirror 
to a pinhole for spatial filtering. The reference wave is 
collinearly recombined with the vortex beam and the re¬ 
sulting interference pattern is recorded as a function of 
delay with a CCD camera. Since the expanded reference 
beam is much larger than the original beam, we assume 
it to be homogeneous across the CCD array. This way, 
an interferogram is obtained for each pixel, containing 
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the spatially resolved phase information between the ref¬ 
erence and the vortex beam. 

The temporal characterization of the reference wave is 
performed using the d-scan technique [25]. A small cen¬ 
tral portion of the beam is selected by an iris with about 
200 fjbm diameter (see Fig. 1). This portion of the beam 
contains sufficient power for second harmonic generation 
(SHG), while being small enough to be assumed spatially 
homogeneous. An off-axis parabolic mirror focuses the 
beam to a thin (20 fim) KDP crystal and a broadband 
polarizer blocks the fundamental NIR light from entering 
the fiber that sends the SHG signal to a spectrometer. 
The SHG spectrum is recorded as function of glass inser¬ 
tion (dispersion) around the optimal compression inser¬ 
tion. From the SHG trace and the measured fundamen¬ 
tal spectrum, the spectral phase (Fig. 2(c)) is retrieved 
using an iterative numeric algorithm [25]. From the re¬ 
trieved phase, the corresponding time domain pulse can 
be obtained at any glass insertion. The shortest pulses 
obtained have a FWHM duration of 5.5 fs (Fig. 2(d)). 
We assume that the measured spectral phase is a good 
approximation to that of the spherical reference pulse 
used in the main interferometer setup, which is used for 
the spatiotemporal reconstruction presented below. 


| U(t) \ exp [ii/;(i)\, or its corresponding complex spectrum 
U{u)) = \U(uj)\exp[i(f)((jj)], interfering with a reference 
pulse U r (t) = \U r (t)\ exp[i'0 r (t)]. The obtained linear 
interferogram, measured with an intensity detector as a 
function of delay, is given by 

I( T )= j \U{t) + U r (t-T)\ 2 dt. (3) 

The Fourier transform of this interferogram can be writ¬ 
ten as: 

ni<?)}=F{j \U(t)\ 2 dt + j \U r (t)\ 2 dt^j 
+ U(uj)U*(u>) 

+ U*{-u)U r (-uj). (4) 

Equation 4 is composed of three terms, the first one 
being the DC component of the interferogram and the 
second and third terms being the complex conjugate of 
each other. The three terms are easily separable in the 
frequency domain. We choose to isolate the second term: 

A(oj) = U(uj)U*(uj) 

= \U(u)\\U r (u)\exp[i^>(cj) — icj) r (uj)]. (5) 
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Fig. 2. Temporal characterization of the reference beam, 
(a) Measured and (b) retrieved d-scan traces, (c) measured 
fundamental spectrum (black line) and retrieved phase (red 
line) and (d) reconstructed temporal profile at the reference 
glass insertion (zero amount of glass on the d-scan plots). 
The pulse duration (FWHM) at the reference insertion is 
5.5 fs FWHM. 


3. Spatially Resolved Fourier Transform Spectrome¬ 
try 

We first explain the principle of Fourier transform 
spectroscopy considering only one point in space, thus 
omitting the spatial coordinates. Consider an ultra- 
short pulse, described by the complex field U(t) = 


If the reference phase is known and the reference 

spectrum \U r (u)\ covers the bandwidth of the field to be 
measured, then U(uj) can be retrieved: 

U(u) = exp[^ r M]. (6) 

\U r (uj)\ 

We now consider the spatial transverse dependence 
of the waves. If the reference wave is homogeneous, 
U r (x,y,uj) = U r (uj), the spatially resolved spectral am¬ 
plitude and field of the unknown pulse can be found: 

U(x,y,u>) = exp[^ r (u)], (7) 

\Ur(u) I 

as well as the corresponding time domain field U(x, y , £), 
by Fourier transformation. In practice, it is much eas¬ 
ier to produce a spherical reference wave (Fig. 1). The 
spherical curvature of the reference wave can easily be 
removed numerically from the measured interferogram. 

It should be noted that the spatial phase front can be 
extracted for every wavelength even without knowledge 
of the spectral phase of the reference. However, the com¬ 
plete spectral phase and thus the pulse duration of the 
unknown pulse cannot be obtained without knowledge 
of the reference pulse. 

4. Frequency-Dependent Spatial Characterization 

From the complete spatiotemporally characterized field 
U (x, y, u)) the frequency-dependent characteristics of the 
vortex can be investigated. For any given wavelength A a 
(or corresponding frequency uj a ) the complex 2D array 
U{x^y^u) a ) can be selected and the intensity (the abso¬ 
lute square of the amplitude) and phase can be obtained. 
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Fig. 3. (a-c) Intensity and (d-f) phase profiles for three dif¬ 
ferent wavelengths. The dimensions of each of the plots are 
4 mm x 4 mm. 


Fig. 3 shows intensity and phase profiles at the CCD 
detector for three different wavelengths. For the design 
wavelength of the phase plate (800 nm) a clean vortex is 
obtained. The effect of the plate’s chromaticity increases 
as the wavelength departs from the nominal one. 

At 631 nm (Fig. 3(a)) strong diffraction effects, giving 
rise to horizontal stripes, are noticeable. They originate 
from the phase jump introduced by the plate (« 2.57T 
at 631 nm). At longer wavelengths the effect is not 
as clearly pronounced, mostly because material disper¬ 
sion is lower making the phase jump closer to 2i r. Fig¬ 
ures 3(c,f) show the reconstructed profile at 952 nm. 
At this wavelength, the phase jump is estimated to be 
around 1.77T, leading to less evident distortion on the 
intensity and phase profiles. We find that for any wave¬ 
length within the spectrum of the pulse, there is a phase 
singularity present in the center of the beam with the 
intensity dropping to zero. 

The effects originating from phase plates with phase 
jumps not equal to 2tt (see e.g. [26]) can be reproduced 
by simple propagation simulations. Figure 4 shows an 
optical field, originally Gaussian, just after a spiral phase 
plate with a phase jump of 2.5tt (Fig. 4(a,d)) and after 
propagation to the position of the detector (Fig. 4(b,e)). 
The simulated intensity and phase profiles (Fig. 4(b,e)) 
qualitatively agree with the measured profiles at 631 nm 
(Fig. 3(a,d)). The differences between measured and 
simulated phase profiles away from the center are mostly 
due to astigmatism which was not accounted for in the 
simulation. 

We also back-propagated the measured field 
(Fig. 3(a,d)) to the plane containing the phase 
plate. The results are plotted in Fig. 4(c,f), clearly 
showing a phase discontinuity (around 0.57r) due to the 
thickness jump of the phase-plate, as well as a slight 
astigmatism of the used beam. 




Fig. 4. Simulated intensity and phase profiles (a,d) after 
the spiral phase plate and (b,e) numerically propagated to 
the position of the detector. We assume a Gaussian beam 
and a phase step of 2.57T due to the plate. (c,f) Calculated 
intensity and phase profiles obtained by back-propagating 
the measured field (Fig. 3(a,d)) to the plane of the plate. 
The dimensions of each of the plots are 4 mm x 4 mm. 


5. Time-Dependent Spatial Characterization 

From the complete description in the spectral domain 
U(x,y,uj) the spatiotemporal field U(x,y,t) can be ob¬ 
tained by simply Fourier transforming at each point in 
space. To achieve the shortest possible vortices as shown 
here, 1 mm of fused silica (corresponding to the sub¬ 
strate of the spiral phase plate) was removed by adjust¬ 
ing the wedges used for compression. 

Figure 5(a) shows the spatiotemporal intensity profile 
\U(x, y, t)\ 2 as an isosurface set at half the peak inten¬ 
sity value. The projections on the vertical and horizon¬ 
tal planes in Fig. 5 are obtained by taking cuts through 
the center of the beam. The time duration of the pulse 
taken at any individual point on the main ring is ap¬ 
proximately the same as for the reference pulse, i.e. « 
5.5 fs. However, considering the vortex as a whole, the 
duration is longer, around 8 fs. The largest contribu¬ 
tion to that is the phase plate itself, which delays parts 
of the beam with respect to others. Furthermore, resid¬ 
ual astigmatism and diffraction effects resulting from the 
chromaticity of the phase plate contribute to the longer 
pulse duration. The phase singularity, apparent as van¬ 
ishing intensity in the center of the beam, can clearly be 
observed (inset of Fig. 5(a)). 

To visualize the spiraling nature of the electric field 
in time we take the real part of the complex field, 
E(x,y,t) = SR{C7(x,y,£)}. Figure 5(b) shows isosur¬ 
faces of E(x,y,t) set at a value of zero, which shows the 
wavefront of the spiraling field, while Fig. 5(c) shows 
isosurfaces set at 0.7 and -0.7 of the peak field value, 
respectively. Note that the absolute phase of the elec¬ 
tric field is arbitrary, since the carrier-to-envelope offset 
phase of the reference was not measured. 
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Fig. 5. Three dimensional reconstruction of the vortex beam, (a) Intensity profile in space and time obtained by taking the 
square of the absolute field. The plot shows an isosurface set at half the peak intensity. The inset shows the same surface 
from a view along the propagation direction, visualizing the phase singularity (dimensions of inset plot are 4 mm x 4 mm), (b) 
Isosurface set at 0 of the real part of the complex field. For clarity, only regions where the intensity is higher than 0.4 times 
the peak intensity are shown, (c) Isosurfaces of the real part of the field, set at 0.7 (red) and -0.7 (blue) of the peak field. 


6. Conclusions 

We have demonstrated the generation and complete 
characterization of ultrashort vortex pulses, using simple 
to implement yet powerful techniques. 

The combination of the d-scan technique with spa¬ 
tially resolved Fourier transform interferometry has 
proven to be a powerful technique for the characteri¬ 
zation of ultrashort optical pulses with complex spatial 
distribution and spatiotemporal couplings, being espe¬ 
cially suited for measurements requiring high spatial res¬ 
olution with accurate phase sensitivity across the whole 
beam profile. 

The spiral phase plate, although not designed for 
broadband operation, has a good trade off between per¬ 
formance and ease of use, since it is extremely simple to 
align, has low losses, and is well suited for high-power ap¬ 
plications. Despite the chromatic response of the spiral 
phase plate, which leads to a break of circular symmetry 
for spectral portions of the pulse away from the design 
wavelength, the phase singularity in the center of the 
beam remains intact. 

Our characterization shows that high power ultrashort 
vortex pulses can conveniently be generated with spiral 
glass phase plates. This method should be useful for 
strong field physics with singular light fields. 
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